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INTRODUCTION 


riu'  ]inrposc  of'  this  report  is  t(>  propose  a microscopic  model  of 
tlurmal  initiation  of  a solid  explosive  and  to  examine  existing  data 
foT-  some  simple  solid  explosives  to  test  the  validity  of  the  proposal. 
I’a.st  el'fT)rts  to  understand  the  [irocess  in  solids  have  involved  macro- 
scopic formulations  using  continuum  mechanics  and  emjihasicing  pro- 
cesses in  the  detonation  regime  (Ref  1 and  21.  However,  these  ap  - 
proachc.s  are  not  fundamental  in  the  sense  that  tliey  do  not  view  the 
solid  as  an  array  of  (piantum  jtarticles.  'Hie  initiation  of  explosion, 
the  jii'ocess  wherehs'  a chemical  reaction  commences  which  then  induces 
a self-sustaining  series  of  exothermic  reactions,  is  operative  in  the 
sol  ill  state  and  therefore  should  be  influenced  by  the  solid  state 
properties  of  the  material. 

It  will  bo  proposed  that  the  inititil  reaction  in  the  explosive 
solid  is  caused  by  direct  contact  of  adjacent  anions  of  the  solid 
which  can  arise  because  the  modes  of  the  explosive  lattice  are  softer 
and  the  potentials  of  the  lattice  more  anharmonic.  Reaction  rates 
for  such  a mechanism  will  be  developed  in  terms  of  the  force  con- 
stants of  the  normal  modes  of  vibration  for  the  case  of  a simple 
ciiliic  lattice.  The  resulting  reaction  rates  will  be  used  to  examine 
the  effect  of  softer  and  more  anharmonic  modes  on  the  reaction  rates. 
Available  measurements  of  lattice  mode  frequencies  in  [lotassium  azide 
(K.N'O,  a non-explosive,  and  the  isostructural  explosive  tha’lous 
azide  (TI.N3)  will  be  compared  to  test  the  validity  of  the  model. 

Other  j)i'operties  of  the  materials  will  also  be  considered. 


THEORY 


Ihe  initiating  reaction  in  the  explosion  of  a solid  is  consider- 
ed to  arise  by  the  chemical  interaction  of  two  or  more  anions  of  the 
solid.  The  explosive  solid  then  must  have  ions  capable  of  'nteract- 
ing  exothermically.  The  reaction  is  not  believed  to  involve  direct 
dissociation  of  the  anion.  .Such  a process  is  endotliermic  and  of  it- 
self could  not  lead  to  detonation.  TTie  products  of  this  decomposi- 
tion could,  however,  react  with  each  other  exothermically.  However, 
at  .TOn°(I  (the  typical  initiating  temperature  of  primary  explosives) 
or  less,  it  can  be  estim.ated  that  there  would  be  an  insufficient 
number  of  j)roducts  available  which  would  lead  to  a reaction  with  suf- 
ficient yield  of  energy  to  cause  detonation.  It  seems  more  probable 
that  the  reaction  is  between  existing  ions  of  the  lattice  rather 
than  products  of  the  direct  decomposition  of  the  anions  of  the  solid. 
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11'  this  IS  tile  case,  the  question  arises  as  to  how  tlie  ions  of  the 
lattice  can  jtot  sufficiently  close  to  undergo  chemical  reaction  at 
the  initiation  temperature.  It  is  proposed  that  the  chemical  contact 
hetween  the  ions  of  the  lattice  is  induced  by  large  amplitude  vibra- 
tions of  the  ions  at  the  initiation  temperature.  There  are  two  as- 
pects of  the  potential  that  influence  the  ions  ability  to  make  con- 
tact with  the  adjacent  ions.  If  an  ion  at  position  r^  sees  a poten- 
tial of  the  form, 


V(r  - r^)  = K(r  - rg)^  - gfr  - Tq) ^ (1) 

the  average  displacement  <r  - tQ>  (not  the  amplitude  of  vibration) 
for  the  case  where  kT  is  large  compared  to  the  anharmonic  term  is 
(Ref  7>) 


■ r - rQ>  = (3g/4K^)kT 


(2) 


Note  that  the  average  displacement  (but  not  the  amplitude  of  vibra- 
tion) is  zero  for  a harmonic  potential.  The  expression  shows  clearly 
that  when  g is  large  and  K is  small,  the  average  displacement  is  large 
Note  also  that  small  K will  mean  large  amplitudes  of  vibration.  Per- 
haps then  one  reason  why  TIN3  is  explosive  and  KN3  is  not  is  that 
the  explosive  lattice  has  softer,  more  anharmonic  potentials,  which 
serve  to  enhance  the  probability  of  chemical  contact  between  the  ions. 
Only  force  constants  of  normal  modes  which  involve  motions  that  allow 
the  ions  of  the  lattice  to  approach  each  other  need  be  softer.  A 
further  desirable  feature  of  the  model  is  that  in  the  explosive  these 
modes  would  be  more  susceptible  to  softening  with  increasing  tempera- 
ture . 

Relationship  Between  Force  Constants  and  Reaction  Rates 

It  is  assumed  when  the  adjacent  ions  in  a lattice  have  a certain 
critical  amplitude  of  vibration,  A^,,  the  ions  will  be  sufficiently 
close  to  chemically  interact.  Referring  to  Figure  1,  consider  two 
adjacent  ions  of  a lattice  having  mass  mj  and  m^,  separated  by  a dis- 
tance, 1.,  and  vibrating  in  a potential  V(x).  Only  ions  of  potential 
energy  V(A^J  or  above  will  be  able  to  react.  Note  that  the  vibration 
of  the  particles  must  involve  motions  which  allow  the  ions  to  approach 
each  other.  Ttiis  limits  the  specific  modes  of  vibration  that  can  be 
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involved  in  the  reaction.  Tlie  fraction  of  the  particles  in  the 
i th  mode  which  have  energy  V(A^)  or  greater  is 


f^  = exp  (-V(A^)/2kT) 


The  total  number  of  particles  in  the  system  having  amplitudes  greater 
than  Aj,  is  nfi  where  n is  the  total  of  all  particles  in  the  system. 
The  particles  on  the  average  will  approach  each  other  t/2  where  t 
is  the  period  of  oscillation.  Thus  the  reaction  rate,  R,  is 


— nexp 


(-V(A^)/2kT) 
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It  is,  of  course,  a serious  over  simplification  to  assume  that  a re- 
action is  turned  on  suddenly  at  some  precise  separation  or  amplitude 
of  vibration  A . More  realistically,  the  reaction  probability  will 
be  a gradual  function  of  the  separation  of  the  reacting  species. 
Equation  4 can  be  multiplied  by  a factor  which  is  a function  of  A 
and  describes  the  probability  of  reaction  once  a given  A is  achieved. 

Example  of  a Simple  Cubic  Lattice 

The  above  discussion  indicates  how  reaction  rates  in  a solid  can 
depend  on  the  force  constants  and  frequencies  of  the  normal  nodes  of 
vibration  of  the  solid.  In  order  to  show  the  role  of  these  nodes  in 
the  reaction  process,  a simple  face-centered  cubic  explosive  solid 
will  be  considered.  Such  a simple  structure  for  an  ex])los-vc  solid 
is  not  unrealistic.  The  high  temperature  phase  of  T1N3,  from  which 
the  material  explodes,  is  known  to  be  a body-centered  cubic  lattice, 
as  is  ammonium  nitrate  (Ref  4 and  5).  For  purposes  of  illuminating 
the  ideas,  a n’jmber  of  simplifications  will  be  made  in  the  analysis. 
First,  only  nearest  neighbor  interactions  in  the  lattice  will  be  con- 
sidered; also,  the  analysis  will  be  made  for  zero  wave  vector.  In 
the  transverse  optical  mode,  depicted  in  Figure  2,  the  ions  oscillate 
parallel  to  Y,  but  the  wave  is  propagating  along  X.  The  harmonic 
force  constants  are  defined  in  Figure  2.  It  can  be  shown  that  for 
K=0,  and  including  only  nearest  neighbor  interactions,  the  frequency 
is  given  (Ref  6)  by 


m,  + m, 

/ • '"»•  (5) 
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A similar  analysis  can  be  made  for  the  other  modes.  However, 
when  the  motions  of  the  accoustical  modes  are  examined,  it  is  seen 
that  they  do  not  allow  the  ions  to  approach  each  other.  In  this 
example,  the  reaction  is  considered  to  be  between  ions  of  mass  mj 
and  m, . If  the  fact  that  the  ions  of  the  lattice  are  really  polar- 
izable ions  is  included  in  tlie  analysis,  the  eftective  force  constant 
for  the  transverse  optical  mode  can  be  shown  to  be  (Ref  6) 

Kef=  (2A  + 4B  - qV(|^  - a))  (6) 

where  q is  the  charge  on  the  ions,  V is  the  volume  of  the  unit  cell, 
and  u is  the  polarizability  of  the  ion.  ’Hie  important  point  is  that 
for  the  t.o  mode,  the  effect  of  including  polarizable  point  charges 
in  the  analysis  is  to  make  the  force  constant  softer,  for  an  ionic 
lattice,  the  last  term  in  Kpff  can  be  .^0%  of  the  total  force  constant. 
Thus  the  polarizability  of  the  ions  plays  an  important  role  in  de- 
termining the  softness  of  a force  constant,  a fact  wliich  is  germane 
to  explosive  solids.  Combining  liquations  4 and  5,  a reaction  rate 
for  the  probability  of  interaction  of  ions  1 and  2 in  the  harmonic 
approximation  for  K^=0  of  a face-centered  cubic  lattice  is  obtained  to 
to  bo 


'1110  influence  of  the  strength  of  the  force  constant  on  the  reaction 
rate  can  be  seen  by  considering  a hypothetical  situation. 

In  figure  ,3,  the  reaction  rate  versus  force  constant  is  plotted 

for  a number  of  different  assumed  values  of  A where  m,  is  taken  as 

c ^ 

a thallous  atom  and  m2  the  azide  ion.  The  result  indicates  that,  de- 
pending on  the  value  of  A^,  the  reaction  rate  will  increase  with  the 
decreasing  force  constant.  Interestingly  the  mode  need  not  soften 
completely  to  zero  for  the  reaction  rate  to  reach  a maximum.  fhe 
maximum  possible  reaction  rate  for  a given  value  of  A^  can  be  obtained 
from  liquation  7 to  be 

nil  "*2  (kT/A^)  e'l/2 

^max  J mimp 
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If  it  is  assiimotl  that  tho  maximiim  reaction  rate  is  the  rate  jiecessary 
for  explosion,  the  above  expression  indicates  a relationship  between 
the  initiation  temperature  and  the  mass  of  the  ions  in  the  lattice. 

When  the  above  analysis  is  made  including  the  anharmonic  terms 
in  the  potential,  the  bimolccular  reaction  rate  can  be  enhanced  with 
increasing  anharmoni ci ty  and  decreasing  force  constant.  Thus  given 
two  isostructural  materials.  which  is  explosive  and  SXj,  which  is 

stable,  the  lattice  potential  should  lie  softer  and  more  anharmonic  in 
the  EXj^.  Ihese  differences  n the  potential  should  be  reflected  in 
a comparison  of  a number  of  basic  properties  of  the  materials,  such 
as  mode  Grimeisen  parameters,  dielectric  constants,  coefficients  of 
thermal  expansion,  and,  of  course,  the  mode  frequencies  themselves. 


EXPERIMENTAL  EVIDENCES  FOR  THE  MODEL 


In  this  section  the  properties  of  some  primary  explosives  solids 
will  be  examined  to  see  if  there  is  any  possible  experimenta  support 
for  the  proposal.  Although  other  materials  will  be  considered,  the 
emphasis  will  be  on  azides  and  fulminates  because  there  is  more  avail- 
able data  for  these  materials  than  for  other  solid  primary  explosives, 
and  because  they  are  structurally  simple.  In  particular,  structural 
isomorphs  such  as  TIN3  and  K.N3,  have  relatively  simple  structures  and 
quite  strikingly  different  stabilities.  Such  isomorphs  allow  a com- 
parison of  fundamental  properties  to  test  the  ideas  of  the  model. 

Equation  S suggests  that  there  should  be  a relationship  between 
the  initiation  temperature  and  the  cation  mass  in  a series  of  explo- 
sives having  the  same  anion.  Tlie  initiation  temperature  versus  ca- 
tion mass  for  a series  of  monovalent  fulminates  is  plotted  in  Figure  4 
and  indicates  a trend  of  the  initiation  temperature  with  cation  mass. 
Interestingly  a similar  plot  for  the  divalent  azides  shows  a trend, 
but  in  the  opposite  sense.  It  should  be  emphasized  that  the  concept 
of  an  initiation  temperature  is  not  well  defined.  The  major  diffi- 
culty is  that  there  is  a dependence  on  the  time  the  material  is  held 
at  a given  temperature  before  it  initiates,  which  is  related  to  the 
geometry  of  the  sample  and  the  thermal  conductivity  as  well  as  to  the 
method  of  measurement.  To  minimize  these  problems,  only  initiation 
temperatures  obtained  by  one  method  with  one  consistent  definition 
were  used  in  the  plot  (Ref  7) . 
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Despite  the  fact  that  materials  of  different  electronic  and 
crystal  structure  are  being  compared  in  Figure  4,  it  is  remarkable 
that  the  trend  exists  at  all.  Hiis  trend  may  be  suggestive  of  the 
role  of  a lattice  dynamical  mechac'sm  in  the  initiation  process. 

Hie  most  direct  way  to  test  the  assumptions  of  the  model  is  to 
examine  measured  frequencies  of  lattice  modes  in  the  structurally  iso- 
morphic explosive  solids.  Both  TIN^  and  KN3  have  [),J^  symmetry  and  a 
body- cente red  tetragonal  unit  cell  which  is  shown  in  Figure  5 (Ref  8). 
In  Table  1 are  listed  the  measured  K=0  mode  frequencies  for  TIN3  and 
KiN3,  obtained  from  Raman  and  infrared  studies  (Ref  8 and  9j . The  data 
indicates  that  tiie  T(EjjJ  and  fCEgi  transl  tional  modes  and  the  R(Eg) 
librational  mode  of  the  azide  ion  have  s- hstant i al ly  lower  frequency 
in  the  explosive  solid.  Ttic  motion  of  'he  individual  ions  associated 
with  the  modes  that  have  lower  frecpiency  in  the  explosive  are  shown 
in  Figure  6.  It  is  clear  that  the  motion  of  the  ions  for  the  TCF.yl 
and  the  R(l;  ) modes  involves  the  approach  of  the  azide  ions  to  each 
other.  It  Ts  necessary  to  look  at  the  data  a little  closer  to  assess 
whether  a reduced  frequency  necessarily  means  a lower  force  constant. 
It  could  be  simply  a result  o"-'  the  heavier  cation  mass  in  the  TIN3. 
Because  the  force  constant  and  the  mode  frequency  are  simply  related, 
the  R(E  ) mode  will  be  considered  as  an  example.  If  mode  coupling 
and  anhtrmoni ci ty  are  neglected,  the  frequency  of  the  librational 
mode  is  given  by 


u) 
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(9) 


The  second  derivative  of  the  potential  in  the  expression  is  the  spring 
constant  of  the  mode.  The  I is  the  moment  of  inertia  of  the  azide  ion, 
which  is  2mA^.  Tlie  calculated  force  constant  for  the  K(Ea)  mode  of 
TIN3  is  .^.0  X 10"^  dyne-cm  compared  to  25.0  x 10"^  dync-cm  for  KN3. 

Thus  the  force  constant  for  this  mode  in  the  explosive  is  about  a fac- 
tor of  eight  softer  in  the  explosive.  Modes  with  softer  force  con- 
stants will  have  larger  amplitudes  of  vibration.  In  Table  2 are 
listed  amplitudes  of  vibration  of  the  Eg  librational  mode  in  TIN3  and 
KN3  at  293  K obtained  from  neutron  diffraction  and  Raman  data  (Ref  8, 

9,  and  10).  It  is  clearly  seen  that  the  N3  in  TIN3  will  have  a larger 
amplitude  of  libration  at  a given  temperature  than  in  KN3.  Of  course 
it  is  the  amplitude  at  the  initiation  temperature  that  is  important. 

A rough  estimate  of  the  amplitude  of  libration  at  600  K can  be  made 
from  Crui ckshank ' s foinrila  for  the  case  where  kT  > hv  (Ref  II) 


0^  = kT/4  TT^Iv^ 


(10) 


6 


I rom  ,1  va i 1 ah  I tcnijicratiirc'-deperulent  moasuromcnts  of  the  fre- 
quency of  the  ino(J(>,  the  frequency  of  the  mode  at  600  K can  he  obtained 
nie  calculation  leads  to  a Lihrational  amplitude  of  2(i  decrees  in  TIN^ 
comjiared  to  8 de.’.rees  in  KN3.  'ITius  there  is  a siftnificant  difference 
in  the  ;implitiules  of  vibration  for  this  mode  between  the  explosive  and 
the  nonexplosive.  Certainly  such  a large  amplitude  of  vibration  would 
allow  adjacent  azides  to  make  a closer  approach  to  each  other  at  oOO  K 
in  TIN,.  'Hie  lihrational  motion  of  the  azide  ion  associated  with  the 
K(l.^)  mode  involves  an  oscillation  of  *he  ion  about  110  (or  ITO  for 
the  other  site).  Such  a motion  would  bring  the  end  of  one  azide  ion 
(which  are  adjticent  to  each  other  along  the  C axis)  closer  to  the 
center  of  the  nearest  neighbor  along  C.  faking  the  covalent  radii  of 
the  nitrogen  as  1.5  A and  using  the  lattice  parameters  it  the  highest 
temperature  of  measurement  (498  KJ,  it  is  calculated  that  an  angl«'  of 
48  degrees  would  allow  the  nitrogens  of  the  azides  to  "touch."  Hie 
force  constant  at  the  initiation  temperature  can  be  calculated  from 
liquation  9,  then  used  in  liquation  7>  to  estimate  the  fraction  of  the 
total  number  of  azides  in  the  lattice  that  would  have  amplitudes  large 
enough  to  allow  contact.  'ITie  result  indicates  that  more  t’lan  2.5°o  of 
the  ions  of  the  lattice  would  have  the  necessai'y  amplitude,  a surpris- 
ingly large  number. 

It  must  be  emphasized  that  there  is  no  evidence  to  conc’ude  that 
such  a contact  of  the  covalent  i-adii  as  described  above  will  induce 
chemical  reaction  between  the  azides.  The  motivation  of  the  discus- 
sion has  been  to  demonstrate  that  it  is  quite  reasonable  that  a sig- 
nificantly large  percentage  of  the  azides  can  get  closer  in  the  e.xplo- 
sive  because  of  the  softer  lihrational  mode.  Certainly  contact  of 
covalent  radii  is  not  a valid  criteria  for  inducing  chemical  reaction. 
A similar  analysis  for  the  translational  mode  leads  to  the  same  con- 
clusion, that  a s'gnificant  number  of  the  ions  of  the  lattice  can  get 
closer  to  each  other  in  the  explosive  compared  to  the  nonexplosive. 
Further  research  is  necessary  to  determine  which  of  the  modes  is  more 
conducive  to  causing  reaction  between  the  azide  ions.  Such  an  inves- 
tigation would  reciuire  a study  of  the  reaction  probability  as  a func- 
tion of  the  different  ways  the  two  ions  can  approach  each  other. 

In  this  discussion  of  the  libratlonal  mode,  the  role  of  anharmo- 
nicity  was  ignored.  If  it  is  included,  the  effect  could  be  to  in- 
crease the  susceptibility  of  the  ions  to  get  closer  as  pi'eviously  dis- 
cussed relative  to  liquation  2.  A measure  of  the  anharmonic  term  in 
a potential  of  a mode  of  vibration  is  obtained  from  the  mode  (Iruneisen 
constant  y j,  which  is  obtained  by  pressure  measurements  of  the  mode 
frequencies  in  conjunction  with  a knowleilge  of  the  coefficient  of 
thernial  expansion,  iiie  V . lias  been  measured  for  the  K(!i,)  mode  in 
KN3  to  be  2.04,  comjiarcd  '^o  6.5  in  ’I'lN,,  clearly  indicatVng  a larger 
anharmonic  contribution  to  tl'.e  a lattice  potential  of  the  explosive 


lattice  (Kef  12).  Ilie  fnineisen  parameters  have  not  been  reported 
for  the  translational  modes  in  these  materials.  Another  indicator  of 
a higher  anharmoni  ci  ty  in  TIN.^  is  the  larger  coefficient  of  thermal 
expansion.  At  a 180°C  the  expansion  coefficient  along  the  C axis  is 
17. .S  x lO"'*  for  TIN3  compared  to  9.2  x 10”'^  for  KN^  (Kef  12). 

With  one  exception,  no  other  solid  explosive  which  can  he  ther- 
mally initiated  from  the  solid  state  has  been  the  subject  of  sufficient 
lattice  dynamical  studies  to  allow  a test  of  the  model.  The  one  ex- 
ception is  a materi;.!  quite  different  from  the  azides,  S^N4.  The 
material  can  be  initiated  at  20S°C  and  pressure  and  temperature  studies 
of  the  mode  frequencies  have  been  made  (Kef  13) . The  substance  has  a 
monoclinic  unit  cell  with  four  S^N,,  molecular  units  per  unit  cell 
(Ref  14)  Tlie  internal  structure  of  the  molecule  is  that  of  a 

cradle-shaped  ring.  The  lowest  observed  mode  at  4.3  cm’’  (at  30  K) 
is  a lihrational  mode  which  softens  by  almost  50%  as  the  initiation 
temperature  is  apjiroached  (Ref  13).  Ihe  mode  is  also  sensitive  to 
pressure.  Ihe  large  amplitude  associated  with  this  libration  would 
induce  contact  between  two  S4.N1,  units.  Kemmey  (Kef  l.S)  has  calculated 
the  potential  energy  of  the  linear  azide  ion  as  a function  of  the  dis- 
placement of  the  ion  along  its  axis  in  l’b(N^)^.  Ihc  potential  )f  the 
azide  ion  at  one  site  was  found  to  be  quite  soft  and  anharmonic. 


CONCLUSIONS 


A mechanism  of  thermal  initiation  of  a solid  explosive  has  been 
considered.  It  is  proposed  that  chemical  reaction  leading  to  detona- 
tion occurs  because  adjacent  anions  at  the  initiation  temperature  can 
make  chemical  contact  as  a result  of  having  large  amplitude  vibrations. 
Large  amplitude  vibrations  would  result  from  softer  modes  of  vibra- 
tion. Only  force  constants  of  modes  which  allow  the  ions  to  approach 
need  be  softer.  Ihe  model  does  explain  how  a large  number  of  reac- 
tions can  occur  in  a relatively  short  time.  Because  of  the  absence 
of  experimental  data,  no  definitive  quantitative  test  of  the  model  is 
possible  at  this  time.  However,  an  examination  of  mode  frequencies 
in  Tl.Vj  and  KNj  appears  to  support  the  possibility.  Modes  which  allow 
the  azides  to  approach  in  TIN3  appear  to  he  softer  in  the  explosive. 

The  proposal  has  been  formulated  in  the  simjilest  possible  terms.  A 
more  vigorous  development  is  not  warranted  at  this  time  due  to  the 
sparcity  of  experimental  data  to  test  the  mechanism. 
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Table  1 


Frequencies  of  lattice  modes  of  KN3  and  TIN, 
at  K 0 measured  at  room  temperature  from  Raman  spectroscopy 


T(*2u> 

R(B,  ) 

<g 

cm‘  * 

R(E  ) 
9 

Ref 

KN3 

1 66 

80 

192 

159 

k47 

8.9 

130 

T1N3 

128 

83 

36 

173 

51 

8,9 

96 


Table  2 


Amplitudes  of  libration  of  the  E mode  in  KN,  and  TIN, 
obtained  from  Raman  spectroscopy  and  neutron  diffraction 


Temperature 

0 

0 

(K) 

(Neutron) 

(Raman ) 

Ref 

KN3 

293 

3.51 

3.15 

10 

TIN, 

293 

9.42 

14.58 

10 

J 1 


papk  nr  ' '*  ' 

V » 1^.  friMHID 

— A-  ^ ^ . 


(seconds)  'x10 


j I I — I — 

12  3 4 

dyne/cm  x 10^ 


Fig  3 Reaction  rates  versus  force  constant  of  transverse 
optical  nxide  in  face- centered  cubic  lattice  for 
different  value  of  "Contact  amplitude"  A^. 
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Fig  4 Explosive  initiation  to  temperature  versus  ratio  of  cation  mass  to  cation 
radii  for  a series  of  monovalent  explosive  fulminates 
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